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ABSTRACT 

Context. Recent lunar crater studies have revealed an asymmetric distribution of rayed craters 
on the lunar surface. The asymmetry is related to the synchronous rotation of the Moon: there is a 
higher density of rayed craters on the leading hemisphere compared with the trailing hemisphere. 
Oh' Rayed craters represent generally the youngest impacts. 

Aims. The purpose of this paper is to test the hypotheses that (i) the population of Near-Earth 

5— i 1 

■j-J , asteroids (NEAs) is the source of the impactors that have made the rayed craters, and (ii) that 

C/3 . 

impacts by this projectile population account quantitatively for the observed asymmetry. 
Methods. We carried out numerical simulations of the orbital evolution of a large number of test 
particles representing NEAs in order to determine directly their impact flux on the Moon. The 
f*"^ ' simulations were done in two stages. In the first stage we obtained encounter statistics of NEAs 

on the Earth's activity sphere. In the second stage we calculated the direct impact flux of the 
■ encountering particles on the surface of the Moon; the latter calculations were confined within 

. the activity sphere of the Earth. A steady-state synthetic population of NEAs was generated from 

a debiased orbital distribution of the known NEAs. 

ON . 

Results. We find that the near-Earth asteroids do have an asymmetry in their impact flux on the 
Moon: apex-to-antapex ratio of 1.32 ± 0.01. However, the observed rayed crater distribution's 
asymmetry is significantly more pronounced: apex-to-antapex ratio of 1.65 ± 0.16. Our results 
5^ ' suggest the existence of an undetected population of slower (low impact velocity) projectiles, 



C3 



such as a population of objects nearly coorbiting with Earth; more observational study of young 
lunar craters is needed to secure this conclusion. 

Key words, asteroids - impacts - craters - Moon 

1. Introduction 

It is well known that many satellites of the solar system planets are locked in synchronous rotation 
— their mean rotational angular speed and mean orbital motion is in a 1:1 commensurability. The 
synchronous rotation of these satellites leads to asymmetric spatial distribution of impact craters 
on these satellites: the leading hemisphere tends to have more craters than the trailing hemisphere. 
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Such leading/trailing asymmetries in crater distrib utions have been observed on the Galilean satel- 
lites of Jupiter and on Nep tune's moon Triton (e.g. 



Shoemaker et al 



1982 



1999; 



Zahnleetal 



Schenk & Sobieszczvk , 



1998). 



Such an asymmetry was recently confirmed on the Moon. A detailed analysis of the Clementine 
750-nm mosai c images has revealed that t here is spatial variation in the density of rayed craters 
on the Moon dMorota & Furumotol 120031) . Lunar rayed craters are morphologically young and 



fresh craters with bri ght rays, generally estimated to be younger than 0.8 billion years old (e.g. 



McEwen et al 



1997). A total of 222 rayed craters larger than 5 km in diameter (D) are identified 



in the study area of about 1.4 x 10 7 km 2 . The average density of rayed craters on the leading side 
of the Moon is found to be substantially higher than that on the trailing side, and the observed ratio 
of crater density ( D > 5 km) at the apex to th at at the antapex is about 1.65. Based on a rough 
analytical estimate. iMorota & Furumotol (12003b conclude that this ratio suggests that recent craters 
on the Moon are formed mainly by the near-Earth asteroids which have lower impact velocities, 
rather than comets that have systematically higher impact velocities. Whether or not these conclu- 
sions are correct, it is true that the ratio of crater densities of the leading side and the trailing side 
of the Moon contains a significant amount of information about the type of projectile populations 
that have created craters on the lunar surface over the past ~ 1 billion years and under what kind of 
dynamical conditions. 

The origin of the leading/trailing asymmetry of impact craters on a synchronously rotating 
planetary satellite is qualitatively explained as follows. Assume that the source of impacts is a 
heliocentric population of small objects on modestly eccentric and inclined orbits. In the rest frame 
of the planet, this population appears almost isotropic and the impact velocity vectors have an 
isotropic distribution. The impact craters asymmetry occurs because the satellite in synchronous 
rotation encounters projectiles more often on its leading side than on its trailing side. Furthermore, 
the average impact velocity of projectiles tends to be larger on the leading side than on the trailing 
side due to the difference of average relative velocities between the projectiles that encounter the 
leading and the trailing sides; this leads to systematically larger craters on the leading hemisphere 
compared with the trailing hemisphere. 

The degree of the asymmetric crater distribution is a function of satellite's orbital velocity and 
the average relative velocity between projectiles and the planet-satellite system. When a satellite 
with synchronous rotation has a large orbital velocity around its mother planet, or when the aver- 
age relative velocity between projectiles and the planet-satellite system is small, the asymmetric 
distribution of craters becomes more pronounced. Smaller orbital velocity of the satellite, or larger 
average relative velocity of projectiles tends to diminish the asymmetry of crater distribution. 

The purpose of the present paper is to quantitatively test the hypothesis that impacts from 
the NEA population (with its currently known properties) account for the observed asymmetric 
crater distribution on the Moon. We do this by carrying out extensive numerical integrations of 
test particles to simulate the impact flux of NEAs. In order to obtain impact statistics and impact 
velocity distribution, we calculate direct impacts of projectiles on the Moon without analytical 
approximation. In Section [2] we describe our dynamical model, our numerical method and our 
choice of initial conditions, and a description of the first stage of our numerical simulations in 
which we trace the dynamical evolution of test particles from their initial locations to the edge of 
Earth's activity sphere; our results on NEA encounters with the Earth's activity sphere are given in 
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section im Section[3]describes the second stage of our numerical simulation in which we trace the 
evolution of particles within the Earth's activity sphere to obtain impact fluxes, impact velocities 
and their spatial distribution on the Moon. Section [4] provides a comparison of our simulation 
results with the observations of the lunar crater record. Section [5] is devoted to discussion of the 
results, and Section [6]to a short sum mary and conclusions . 



While our work was in progress, 



Gallant et al 



(2009) published a study with quite a similar 



motivation to ours, which also yielded a similar numerical result about the expected lunar cra- 
tering asymmetry from NEA impacts. Although a large part of our results overlap, the numerical 
models are different. As we describe in section 2, our numerical model is straightforward and in- 
cludes the orbit evol ution of NEA-like particles from their current orbits up to their impacts, while 
(|2009j)'s study use s the NEA orbits wi thout dynamical evolution. In this regard our 

(2009). We also consider some additional important 



Gallant et al 



Gallant et al 



paper serves as a complement to 
implications of the results that were not discussed previously. 



2. Numerical model 



Our numerical model comprises of two stages. In the first stage, our numerical integrations include 
the eight major planets and the Sun, and a large steady-state number of test particles with NEA- 
like orbits (Fig.Q]). We numerically integrate their orbital evolution for up to 100 million years. 
Throughout these integrations, we record all close encounters of the particles that reach the Earth's 
activity sphere. (Note that in the first stage of calculation the Moon is not included.) We use this 
record in our second stage of numerical simulation, in which we adopt the restricted A^-body model 
consisting of the Earth, the Moon, and the Sun, and cloned test particles within the Earth's activity 
sphere (as described in detail in Section [3]). In the second stage, we do not include the effects of 
any planets save the Earth but we include the Moon's gravity. 

For our first stage numerical simulation we generated a synthetic population of particles with 
orbital elements sim ilar to the "debiased" distributions of near-Earth objects (NEOs) described in 
(120021) . (Note that NEOs are largely composed of NEAs, so we will keep using the 



Bottke et al 



term NEAs rather than NEOs in this paper.) Specifically, we generated 18,000 particl e initial con- 
ditions whose distributions of a, e, and / obey the histograms shown in Figure 12 of 



Bottke et al 



d2002l) which gives the debiased orbital distribution of the near-Earth asteroids of absolute mag- 
nitude H < 18.Q The orbital elements of our synthetic NEA population are shown in Fig.Q] This 
population represents a good snapshot of current orbital distribution of NEAs. Studies of impact 
craters in the inner solar system indicate that there has been a relatively constant supply of im- 
pactors over the past three billio n years which has kep t the impactor flux a round the Earth-Moon 



system close to a stationary state (iMcEwen et al 



population is dominated by NEAs dStrom et al 



1997 



2005). 



Ivanov et al 



20021) . and that this impactor 



For the numerical integration scheme we used the regularized mixed-variable symplectic 
method (Levison & Duncanl 1 19941) . The basic framework of our first stage simulation follows 



Ito & Malhotra 



(2006|). When a test particle approaches within the physical radius of the Sun or 



that of planets, we consider the particle to have collided with that body and lost from the NEA pop- 



1 Note that at this stage of our calculation we do not consider at all the size-frequency distribution (or 
absolute magnitude distribution) of the particles. 
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1.0 




semimajor axis [AU] 

Fig. 1. Initial osculating orbital elements of the NEA-like particles in our numerical model, (a) 
Semimajor axis vs. eccentricity, (b) Semimajor axis vs. orbital inclination. (Inclinations are referred 
to the ecliptic at J2000.0.) 



ulation. Also, when the heliocentric distance of a test particle exceeds 100 AU, the particle is con 
sidered lost. Over the 100 Myr length of the simulation, a large fraction (~ 90%, e.g 



Ito & Malhotra 



(120061) ) of the synthetic population would be expected to be removed in this way, and if this loss 
were not compensated, we would not be able to mimic a steady-state NEA flux. We realize the 
steady-state NEA flux in our numerical simulation as follows: for each "lost" particle, we imme- 
diately introduce in our simulation another particle with the original position and velocity of that 
"lost" particle. This procedure achieves a steady-state population of NEAs in our simulation. In 
particular, we verified that the distribution of lunar impact velocity remains steady throughout the 
simulation timespan. 

2.1. Particle encounters with Earth's activity sphere 

Over the 100 myr simulation of a steady-state swarm of 18,000 particles, we found 3,998 collisions 
with Earth. We note that, although the number of planetary collisions is not large in our first stage 
numerical integrations, there are many more encounters at the planetary activity sphere of Earth. 
The activity sphere, also known as the sphere of influence, has a radius of (m/M) 2 ^ 5 d where m is 
the mass of a planet, M is the mass of the Sun, and d is the heliocentric distance of the planet 



(Danby 



1992). Earth's activity sphere, hereafter denoted r\, is about 144 times as large as the 
Earth's radius. In our first stage numerical simulation, we recorded the encounters of particles at 
the Earth's activity sphere over the 100 Myr integration, and found 42,099,969 encounters. This 
number is large enough to establish a time-dependent orbital distribution function of the particles, 
F(a, e, I, a), CI, l\ f) that can be used to create "clones" of particles in order to increase the reliability 
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Fig. 2. Statistics of NEA-like particle encounters at the Earth's activity sphere, r\. (a) Distribution 
of encounter velocity components {v x , v y , v z ); The horizontal unit is km/s. (b) Surface density of 
encounters at the surface of the Earth's activity sphere as a function of the angle from the leading 
points of each of the x, y, and z directions. The vertical unit is 10 6 km 2 . The inset in the panel 
b schematically shows the coordinate system (x,y,z) adopted for this figure: The Sun always lies 
along the -x direction, +z is the normal to the Earth's orbit, the Earth practically goes toward +y 
direction, and r\ is the radius of Earth's activity sphere. Coordinate of the leading point for the y 
data in the panel is (0, r\, 0). 



of the collision statistics between the particles and the Earth or the Moon, as we describe in the 
next section. 

In Fig.[2]we show time-integrated distribution (over the 100 Myr duration of our simulation) of 
encounter velocity components and encounter density (number of encounters per unit surface area) 
at the Earth's activity sphere of the particles. In our simulation the average encounter velocity of 
the particles at the Earth's activity sphere is 22.5 km/s. We show these distributions along all three 
spatial directions, x, y and z. There are several noteworthy features in these distributions. While 
the encounter velocity distributions with respect to the x and z directions are symmetric, the y- 
direction distributions are markedly asymmetric. Because of its very small orbital eccentricity, the 
Earth's orbital motion is practically along the +y direction. We see the expected effect that more 
particles encounter the Earth's activity sphere on its leading side (from the positive y direction) 
than its trailing side: the fraction of encountered particles having negative v y is larger than that of 
the particles having positive v y (the panel (a)). Consistent with this, we see in the panel (b) that the 
encounter density is higher on the leading hemisphere (0 - 90° angle with respect to the leading 
point in the y direction), and lower on the trailing hemisphere. 
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In the panel (b), we also notice a pattern in the encounter density distributions along the x 
direction: the encounter density of the particles is somewhat higher around the angle ~ 90° (with 
respect to the leading point of the x direction, 0, 0)). The pattern of distribution in the x-direction 
is understood by the following consideration. Because of its very small orbital eccentricity, the 
Earth's orbital motion is practically along the +y direction. Along the y axis, the angle with respect 
to the x- leading point, (ri,0, 0), is 90°; this defines the solar terminator at the Earth's activity 
sphere. Particles that encounter Earth in this region have very small velocity relative to Earth, 
particularly when their random orbital velocity is low; the lower average encounter velocity leads 
to the lower encounter frequency of particles in this region. This accounts for the dip near angle 
90° for the x curves in the panel (a). Moreover, this effect would get smaller when the particle 
population has larger random orbital velocity. 

We also notice a pattern in the encounter density distributions along the z direction in Fig. EJb). 
However, the dips near 90° are of noticeably smaller magnitude in this direction than in the x di- 
rection, and the sign is the opposite: we see a concentration of encounters around 90°. We note that 
the vertical scale height of the NEA-like particle population, approximately given by the average 
(a tan/), is much larger than the radius of the Earth's activity sphere (r\ ~ 0.006 AU). This is why 
the encounter frequency of the particles does not vary as much along the z direction compared with 
the x direction. 

From the above results on the distribution of particle encounters on the Earth's activity sphere 
shown in Fig. [2] it is clear that the Earth receives more NEA impacts on the leading (positive y) 
hemisphere than on the trailing hemi sphere. This asymme try leads to AM/PM asymmetries in NEA 



impact events, which is discussed in lGallant et alj d2009l) in detail 



3. Asymmetric impacts on the Moon 

Using the particle encounters at Earth's activity sphere, we generated cloned particles by per- 
turbing the encounter position r and velocity v of each of the original particles so that their 
initial trajectories at the activity sphere become slightly different: r c i one = (1 + <5 r )f original and 
Vcione = (1 + £v)Voriginai> where 5 r and S v are random numbers in the range [-0.1,0.1]. This pro- 
cedure produces a large number of particles that obey nearly the same orbital distribution function 
as the original particles ("F" in the previous description, see Fig. [2]) but with somewhat different 
paths toward the Earth (and the Moon). 

We repeated this cloning procedure five hundred times from the result of the first stage nu- 
merical integrations, generating 21.049895 billion particle initial conditions on the Earth's activity 
sphere. Using these sets of cloned particles, we performed a second set of numerical integrations, 
this time with the restricted A^-body problem including the Sun, the Earth, the Moon, and the cloned 
test particles. Here we did not include the effect of other planets than the Earth, but we included the 
Moon's gravity. All the cloned particles started near the Earth's activity sphere, and were integrated 
until they hit the Earth or the Moon or went out of the sphere. We used the present orbital elements 
of the Moon with true anomaly randomly chosen from to 360° for each of the 500 sets of clones. 
We employed the regularized mixed-variable symplectic method again with a stepsize of 84.375 
seconds (= 2~ 10 days). (This small step size was arrived at by a process of trial to ensure that even 
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Fig. 3. Statistics of impactors on the Moon and on the Earth, (a) Distribution of impact velocity 
on the Moon (solid line) and on the Earth (dashed line) of the clones, (b) Distribution of impact 
angle on the Moon on the entire surface (solid line; denoted as "Total"), on the leading hemi- 
sphere (dashed line; denoted as "Leading"), and on the trailing hemisphere (dotted line; denoted as 
"Trailing"). 



high velocity particle collisions with the Moon were not missed.) Calculations were carried out in 
the geocentric frame. 

The second stage numerical integrations yielded 1,509,364 collisions with the Earth and 73,923 
collisions with the Moon. Fig.[3]shows the distribution of impact velocities and impact angles on the 
Earth and on the Moon. Overall, the average impact velocities of the clones on the lunar surface, 
(~ 22.4 km/s) is almost the same as the average encounter velocity of the original particles at 
the Earth's activity sphere. This means that lunar gravity plays only a minor role in accelerating 
particles to the lunar surface in our numerical model. Not only lunar gravity but the Earth's gravity 
also plays only a small role: average impact velocity of the clones at the Earth's surface is ~ 23.1 
km/s, not being very different from the average impact velocity with the lunar surface, in spite of 
the large difference of the escape velocities from the two bodies (~ 11.2 km/s on the Earth and 
~ 2.4 km/s on the Moon). The ratio of the number of collision s with the Earth and t hose with 
the Moon is found to be 20.4 + 0.1. For comparison, we note that lZahnle & Sleep! (119971) reported 
the ratio of collisional cross sections of the Earth and the Moon as ~ 23, by assuming isotropic 
collisions and average impact velocity of Earth-crossing asteroids to be 16.1 km/s on the Earth. 

Regarding the impact angles on the lunar surface, we note that from simple geometrical consid- 
erations for an isotropic distribution of impact direction, the impact angle distr ibution is expected 
to have a probability density function proportional to sin 2i dShoemaken, 1 19621) where / is impact 
angle and i — » means oblique impact. The results of our simulation, taken all together, are con- 
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Fig. 4. Heliocentric orbital distribution of (a) Initial semimajor axis distribution of all particles 
included in the first stage integrations (solid line), initial semimajor axis distribution of the particles 
(their equivalent clones) that eventually collided with the Moon (dotted line), and final semimajor 
axis distribution of the particles (their equivalent clones) that eventually collided with the Moon 
(dashed line; i.e. when they hit the lunar surface). Note that the dip arou nd 2 AU in the solid line 
distribution is originated from the debiased semimajor axis distribution of J 



Bpttkeetal 



(I2002I) . (b) 



Same as (a), but for eccentricity, (c) Same as (a) and (b), but for orbital inclination. 



sistent with isotropic impact angles (Fig.[3]5). We can mention in passing that there is a small but 
noticeable difference in the impact angle distribution on the leading and trailing hemispheres of 
the Moon: the trailing hemisphere slightly disfavors oblique impacts whereas oblique impacts are 
slightly enhanced on the leading hemisphere, Fig. [3b). This is not of significance for the statistical 
results in the present paper, but it may be of interest for future studies of individual lunar craters. 

It is also interesting to examine the orbital element distribution of the lunar colliders. In Fig.|4j 
we plot histograms of the distribution of a, e, I for the lunar colliders as well as for our synthetic 
NEA initial conditions; for the lunar colliders, we plot histograms of both their initial orbits and 
their orbits just before impact with the Moon. Comparison of the two initial orbit distributions 
shows that the lunar collider population has a higher fraction of low inclinations and low semimajor 
axes compared to the overall NEA initial orbit distribution. Comparison of the initial orbits and final 
(just before impact) orbits of the lunar colliders shows that there is significant dynamical evolution 
of orbital elements that occurs prior to lunar impact: on average, semimajor axes evolve to lower 
values, eccentricities and inclinations evolve to higher values. This evolution takes place during the 
first several million years of their trajectories. 
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The second stage of our numerical simulation yields the spatial distribution of NEA impacts on 
the lunar surface. As mentioned in the introductory section, in order to compare the distribution of 
impacts in our numerical model with the actual lunar crater record, we have to consider a correction 
to the raw numerical results due to the systematic difference in the impact velocities on the leading 
and trailing hemispheres, a difference that owes to the orbital motion of the satellite about its mother 
planet. This correction, which turns out to be quite small for the Moon, arises as follows. 

For a satellite with synchronous rotation, the average impact velocity of projectiles is somewhat 
larger on the leading side than on the trailing side. This difference means that, on average, the 
apparent crater size would be larger on the leading side than on the trailing side (assuming the 
projectile size-frequency distribution (SFD) is not different on the two sides). To illustrate the 
effect this has on the crater densities, consider a power law SFD of craters, as in the solid line 
shown schematically in Fig. [5] where N is the cumulative number of craters per unit area. As a 
result of the higher (lower) average impact velocity on the leading (trailing) side, the SFD curve 
of the craters on the leading (trailing) side becomes shifted toward the positive (negative) direction 
along the horizontal (D) axis, as indicated by the arrow (1) in the figure. This horizontal shift 
is practically equivalent to a vertical shift of the SFD curve toward the dotted line in Fig. [5] as 
indicated by the arrow (2) in the figure, illustrating that N gets larger on the leading side (smaller 
on the trailing side) for the entire range of crater diameter, D. 

The magnitude of the shift depends upon the relationshi p between the impact velocity v 
and the crater size, D. Her e we employ the Pi-group scaling ( iSchmidt & Housen , 



1987 



Melosh. 



1989; 



Housen et al 



1991J) where approximately D oc v® mp with a 



0.44. For the cumulative 
2 which represents well the SFD of the 



SFD of craters, we adopQ N(> D) oc LP with P 
young rayed craters dMorota & Furumotoi 120031) . When the average impact velocity is changed 
from \'o to vi, the cumulative number density of craters at any given diameter D changes from No 
to Ni = No(v\/voyP. This holds for any values of D as long as the crater SFD is expressed by a 
single power law. 

From the results of our second stage simulation, we computed the average impact velocity, 
(vimp) in km/s, of NEAs on the lunar surface as a function of angle from apex, y (degrees), by a 
least squares fit. We find (v; mp ) = -0.00672y+22.7. This indicates that difference of \Vimp) between 
the y = 90° point and the apex (y = 0) or antapex (y = 180°) is less than 0.61 km/s. Compared 
with the average of Vj mp over the entire range of < y < 180°, these velocity differences amount to 
< 2.71%, thus corresponding corrections to the cumulative SFD are given by N]/No ~ 1.02. This 
small difference is owed to the fact that the lunar orbital velocity of ~ 1 km/s is much lower than the 
average impact velocity (v; mp ) of ~ 22 km/s. As a result of this small dependence, apparent change 
of the crater SFD due to the impact velocity difference between the leading side and the trailing side 
is quite modest. This effect, however, would be important when considering the asymmetric crater 
distribution on a satellite with higher orbital velocity aroun d its mother planet such as Ganymede 



around Jupiter dShoemaker et al 



1982; 



Zahnle et al 



2001). 



2 This SFD is also consistent with t he young crater pop ulations linked to impacts by the NEA population 
on all terrestrial planets and the Moon dStrom et alll2005f) . 
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N(>D) 




Fig. 5. Schematic illustration of the change of crater SFD due to apparent change of crater size. D 
is crater diameter, and N(> D) is the cumulative number of craters whose diameter is greater than 
D. Adopted and modified from llshizaki & Furumotol ( ll997l) . 



Including this correction to our second stage simulation, we computed the simulated spatial 
density of NEA impacts on the Moon. Normalizing to unity at antapex, our simulation result for 
the crater density as a function of apex angle are shown in Fig.|6l panel (a). In Fig. [6] we used a 
simple sinusoid with the function form of A + B cos y for a fitting curve where A and B are fitting 
parameters, normalizing A + Bcos 180° = 1. 

For comp arison, panel (b) shows the distribution found from the analysis of observed lunar 
rayed craters dMorota & Furumoto , 2003 ). Note that the n umber of the lunar rayed craters in the ob- 
servational data analyzed by (iMorota & Furumotol 120031) is only 222, while we have about 74,000 
impacts in our simulation. This difference is reflected in the difference of the errorbar magnitudes 
in Fig.|6l which are based on Poisson statistics. 

Examining Fig. [6] what we notice first is that the apex/antapex asymmetry is less prominent in 
the numerical results (panel (a)) compared with the observed lunar rayed crater record (panel (b)). 
The maximum crater density at apex is about 1.65 (normalized to unity at antapex, and estimated 
from the best-fit sinusoid) in the observed crater record, whereas in our simulations, it is 1 .32+0.01 . 

We also examined our numerical model result for trends in the NEA impact density (impact 
flux) with respect to the angles along the x and z axes (Fig. |7). Here we again adopt Poisson 
statistics to assign uncertainties in the data plotted in Fig. [7] We find a tiny dip around the angle 
90° along the x direction (Fig. [7Ja)), although the numerical noise is large. We can interpret this 
dip as related to the dip that we see in the encounter statistics of particles at Earth's activity sphere 
along the x axis (cf. Fig.|2|b)) and as owing to the same dynamical reason. 

Along the z axis, we notice a lower impact density at the polar regions (Fig.|7|b)). We interpret 
this pattern as related to the encounter density at the Earth's activity sphere along the z direction, 
Fig. |2|b). That the number of particle encounters at Earth's activity sphere becomes the smallest 
around the angle = 90° is reflected in the trend found in Fig.|7Jb). The difference in cratering rate 
between at the polar and the equatorial r egions is ~ 10%; this is consistent with the a nalytical esti- 
mate b y 

J2009J). 



Le Feuvre & Wieczorek 



(2008), as well as the numerical result presented in 



Gallant et al 
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Fig. 6. Modeled and observed impact crater distribution on the Moon. We normalized the crater 
density to unity at antapex (y = 180°) using the best fit sinusoid (solid line curve), (a) Numerical 
result including the correction due to the difference of average impact velocity as a function of the 
angular distance (y) from ap ex (y = 0). (b) The observed rayed crater distribution of D > 5 km 



dMorota & Furumoto 



20031) . The inset in (a) illustrates the coordinate system in this frame: The 



Earth always lies along — x direction, the Moon velocity is toward +y direction, and +z is the north 
of the Earth-Moon system. Apex point is defined as (x,y,z) = (0,/?m,0) where Rm is the lunar 
radius. 




60 90 120 

angle from - = R M [deg] 

Fig. 7. Relative impact density (impact flux) on the Moon in our numerical results (a) along x-axis 
with respect to the point (/? M ,0,0), and (b) along z-axis with respect to the point (0,0, Rm)- The 
average impact density over the range of angles to 180° (from x = Rm for the x direction and 
from z = Rm for the z direction) is normalized to unity. Note that the normalization here is different 
from that in Fig. [6] 



5. Discussion 

Does the dynamical model of NEA impacts account for the observed asymmetry of lunar rayed 
craters? The observed crater record has relatively large errorbars compared to our dynamical model 
results (Fig. [6]), and the apex/antapex contrast of the observed and modeled crater densities may be 
considered at least qualitatively consistent with the observations. Such a conclusion would imply 
that the young lunar craters are owed to impacts of the NEAs whose orbital distribution we already 
know. 
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However, the intriguing systematic difference between the degree of asymmetry between our 
dynamical model and the observed crater record, though not enormous, is worthy of comment. 
If we use the best-fit sinusoids in Fig. [6] we can say that the dynamical model accounts for only 
about 49% (= (1.32 - 1)/(1 .65 - 1)) of the observed lunar rayed crater asymmetry. We would like 
to pursue the reasons for this potential discrepancy, hoping to understand and constrain better the 
dynamical origin of the projectiles that have created the lunar rayed craters over the past ~one 
billion years. Because the uncertainties in our numerical model are much smaller than those in the 
observational data (compare the errorbars in Fig. [6^ and Fig.[6j)), we must be very cautious in the 
discussion of this point. We emphasize that the possible explanations discussed below must await 
assessment with improved observational data. 

One possible explanation is related to the impact velocity distributions. The leading/trailing 
asymmetry becomes more prominent when the average relative velocity between the Moon and the 
projectiles is low. The NEA-like particles are, by their dynamical definition, the "slowest" (relative 
to Earth) among all the known small body populations in the solar system. That even these slow 
particles may not fully account for the observed asymmetric distribution in the lunar crater record 
suggests that there may exist a presently-unobserved population of small objects near the Earth's 
orbit that have even lower average relative velocit y than the currently known near-Earth asteroids 



Bottke et al 



(2000 



2002), would not enhance 



do. Conventional debiasing techniques, such as in 
the NEA-like particles with low relative velocity. Rather, such debiasing generally compensates for 
the existence of more distant objects with higher relative velocity (i.e., larger e,I). 

This argument predicts the existence of a hitherto unseen population of slow objects whose he- 
liocentric orbits are close to the Earth-Moon system. We make a rough estimate of the unseen popu- 
lation as follows. The best-fit sinusoid to the observational lunar crater asymm etry is consistent with 
an im pactor population with average lunar impact velocity of 10-12 km/s dMorota & Furumoto , 



2003). Consider an impactor population having a similar shape of the impact velocity distribution 
function as the simulated NEAs (Fig. [3^) but with (vi mp ) = 10-12 km/s rather than the ~22 km/s 
that we found in our dynamical model based on the known NEAs. In such a population, the frac- 
tion of slow objects would be roughly 50% more than the fraction of slow objects in the currently 
known population of NEAs; here we define "slow" NEAs as those having potential lunar impact 
velocity < 12 km/s; such objects would be nearly coorbiting with Earth. In other words, our rough 
estimate of the actual slow NEA population is ~ 50% more than the known slow NEAs. 

Such a population could have remained undetected in observational surveys to date either be- 
cause the surveys have low sensitivity to their orbital parameters or because these objects are fainter 
(smaller and/or darker). The rayed crater record in Fig. [6jb) contains craters with diameter D > 5 
km. On the lunar surface, a crater with D a 5 km can be created by an asteroidal projectile with 
^projectile = 0.2-0.3 km even when the impact velocity is as small as 10 km/s and when the projec- 
tile density is that of porous rock (~ 1.5 g/cm 3 ). These small and slow objects, if they exist in the 
greater numbers that our study indicates, could account for the discrepancy between our numerical 
result and the observed asymmetric crater distribution. More complete observational surveys of the 
near-Earth asteroids can test our prediction. Future progress in the reconstruction of the true orbital 
distribution of NEAs by debiasing techniques would also be useful. 

However, a challenge with the above explanation is that it is not easy to keep NEAs' relative 
velocity too low. There are many complicated resonances in the orbital zones of the terrestrial 
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1997). Thus 



planets that can pump up the random velocities of small bodies (Michel & Froeschld, , 
the slow population would need to be immune to these excitation mechanisms or to be continuously 
resupplied. 

A different explanation might be that larger NEAs with very low relative velocity get frag- 
mented due to Earth's tidal force when they approach the Earth-Moon system. This process would 
increase the number of projectiles (and reduce their size), and may contribute in enhancing the 
asymmetric distribution of craters if the fragments keep the low relative velocity of the parent body 
until they collide with the Moon. 

A third possibility is that the lunar orbital velocity has been larger in the past. A billion years ag o 



the lunar semimajor axis may have been as small as ~ 90% of the current value dBills et al 



1999), 



and has gradually increased to the current value due to the tidal interaction with the Earth. When 
the lunar semimajor axis was 90% of the current value, the lunar orbital velocity with respect to the 
Earth was 17% larger than the current value. This enhancement of the lunar orbital velocity could 
enhance the asymmetric distribution of impacts. But the magnitude of this effect would be limited. 
Even if the lunar orbital velocity has been 17% larger throughout the past one billion years, it would 
be only as small as 1.2 km/s, still too small compared to the average impact velocity. This larger 
value would still be insufficient to explain the difference of the degree of asymmetric distribution 
of the actual lunar craters record (apex/antapex ratio ~ 1.65) from that of our numerical result 
(~ 1.32) when w e consider the semi-analytic estimate of the crater ing rate as a function of the lunar 

hoo% 



orbital velocity dZahnle et al 



2001 



Morota & Furumoto 



20031) . Similarly, 



Gallant et al 



performed a series of numerical simulations in order to check the effect of the smaller Earth-Moon 
distance, and found only a tiny change in hemispherical crater ratio for Earth-Moon distance as low 
as ~90% of the current value (~54 Earth radii). We must note, however, that the history of the lunar 
orbit is thus far predominantly based on theoretical models and is not especially well constrained 
by observations; there may exist an exciting possibility to place an observational constraint on the 

lunar orbital evolution by detailed modeling of the asymmetric lunar crater record. 

A fourth possibility is that the observational lunar crater data of iMorota & Furumotd (12003b 



is incomplete. This crater data consists only of 222 craters of diameter D > 5 km covering about 
a third of the entire lunar surface. There is certainly room for improvement of this dataset. A 
recent brief repor t of a more extensive s earch for lunar rayed craters as small as 0.5 km diameter 



is not conclusive (iWerner & M edvedev 



20101) . A potentially important source of uncertainty and 



confusion in interpreting the spatial patterns in the lunar rayed craters is the uncertain ages of 
these craters. For example, some of the larger craters are argu ed to be older tha n what they had 



been thought fro m a study of optical maturity of their ejecta dGrier et al 



2001). Similar issues 



are pointed out by I Werner & Medvedev! (120101) . Thus, it is possible that the discrepancy could be 
removed with a future complete and correct da taset of young lunar craters. 



Finally, we should comment on the study of 



Gallant et al 



(120091) which already reported results 
that are ove rall rather similar to ours. The first and the largest difference of our numerical model 
from that of ( 



Gallant et al 



( 2009) is that the latter simulated the lunar impacts of NEAs from a sy n- 

(boOl's 



Bottke et al 



thetic NEA sample having fixed orbital elements with distribution following[ 
debiased NEA population (referred to as "source orbits", restricted to 16,307 in the Earth-crossing 
region), whereas we directly integrated the orbital evolution of Bottke's debiased NEA population 
particles (in our first stage simulation) with a steady-state dynamical model. Then, using symmetric 
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characteristics of the orbits, 



Gallant et al 



(2009) effectively multiplied the number of source orbits 
four times. Next, for each of the source orbits, a disk of 10 5 particles with an identical initial ve- 
locity toward the Earth-Moon syst em was created, and integrated until the disk particles reach the 
Earth. Schematic figure of Fig. 4 in[( 



Gallant et al 



(12009b explains well their numerical model of the 
disk; the total number of particles in their numerical model is more than 1 .2 x 10 1 1 . In comparison, 
our model contains a smaller number of particles (0.21 x 10" clones), but we have included plan- 
eta ry perturbatio n s and o rbital evolution in direct numerical integrations of the NEA-like particles. 



Gallant et al 



(2009) concluded that their numerical result of the apex/antapex asymmetry 
(1.28 + 0.01 when considering the ratio of craters within 30° of the apex to those within 30° of 
the ant apex) is completely consiste nt with the value of about 1.6 + 0.1 found in the available crater 
data of iMorota & Furumotd (12003b : they attributed the difference to the large uncertainties in the 



crater data. As we noted above, the differen ce between the dyna mical model and the observations 



is not huge, but, somewhat differently than 



Gallant et al 



(120091) . we conclude that the the results 



of our numerical simulation are only marginally consistent with the observed asymmetry, and we 
ha ve therefore di s cussed at some length several explanations for the possible discrepancy. 

Gallant et al.l d2009l) also pointed out that the impact velocity of the cratering projectiles is 
~ 20 km/s, somewhat higher than values commonly adopted in previous studies, and that this 
has ramifications for pro posed matches betw een the lunar crater size-frequency distributions and 
asteroidal impactors (e.g. lStrom et al.U2005l) . Our calc ulations find the aver age impact velocity to 

d2009l) . Not only the average 



Gallant et al 



be 22.4 km/s, which is even slightly higher than that of 
impact velocity, but also the shape o f the impact velocity distribution in our model is noticeably 
different than in 



Gallant et al 



(200_9|): it is more symmetric about the mean value in our case. We 
attribute these differences partly to the statistical variations of the initial conditions (both studies 
used Bottke et al. (2002)'s debiased a, e, I distribution of NEO orbits, but the particular realizations 
of the set of initial conditions were done independen tly), and partly to the evolution of the particles 
prior to impact on the Moon (IGallant et al.l d2009l) 's numerical model does not account for this 



orbital evolution, but our dynamical model includes this effect). Note that even though the average 
impact veloc ity as well as the sh ape of the impact velocity distribution in our model are different 



from those in 



Gallant et al 



( 2009), the resulting asymmetry in the numerical lunar cratering is quite 
similar to each other. 

Regarding implications for origins of crater populations, certainly the higher average impact ve- 
locity of NEAs on the Moon calls for an update of such studies, but with the caveat that this higher 
value is derived from the currently recognized "debiased" NEA population (which may potentially 
be missing a significant fraction of slow N EAs near the Earth- Moon system). We also note that the 
Late Heavy Bombardment projectiles that , 



Strom et al 



(2005) proposed were not NEAs with the 
steady-state flux but main belt asteroids directly transported from main belt resonance zones to the 
inner sol ar system; such i mpact ors would have an average impact velocity on the M oon of about 
18 km/s Jlto & Malhotra , 



2006), quite similar to the value of 17 km/s adopted in the 



Strom et al 



(200.5) study. 



Currently several lunar missions are underway by several countries (e.g. iNormille & Baglal 



2007). They will yield improved datasets to better determine the nature of the asymmetric distri- 



bution of young craters on the Moon. On the theoretical side, it would be important to improve 
the dynamical models by including more complete physics (such as non-gravitational forces that 
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may be significant in the orbital evolution of small NEAs), and to improve the model estimates of 
observational biases in the known NEA population particularly for those with orbital parameters 
similar to Earth and that are difficult to observe due to their low solar elongation angles. 



6. Summary and Conclusions 

We simulated numerically the spatial distribution of impacts of near-Earth objects, using a nu- 
merical model with a steady-state population of impactors based on current estimates of t he deb i- 
ased near-Earth objects' orbital distribution (as provided by the model of 



Bottkeetal 



(2002)). 



We compared the results of the simulation with the observed asymmetry of the population of 
rayed craters on the leading/trailing hemispheres of the Moon (as measured by 



Morota & Furumoto 



(2003)). Our results and conclusions are summarized as follows. 



1. Our numerical simulation yields a leading/trailing hemisph erical ratio of 1.32 + 0.01 for lunar 



impacts by near-Earth objects. This result is similar to that of lGallant et alj ( 120091) who obtained 
the value 1.28 + 0.01 from a different numerical model. 
2. This result of our numeri cal simulation is only mar ginally compatible with the observed ratio 



of 1.65 ± 0.16 found by Mor ota & Fur umoto (2003). A possible explanation for the discrep 



ancy is that there exists a hitherto undetected population of small objects in heliocentric orbits 
nearly coorbiting with Earth, whose average impact velocities on the Moon are much lower 
than the average impact velocity of the known near-Earth object population. Other explana- 
tions are possible, including the possibility that a more comprehensive study of young lunar 
craters could reveal a smaller leading/trailing asymmetry and thereby remove the discrepancy 
with the dynamical modeling. 

3. The average impact velocity of near-Earth objects on the Moon is found to be 22.4 km/s; Fig. [3] 
plots the impact velocity distribution. 

4. Overall, the impact angles are isotropically distributed, but there is a noticeable deficit of 
oblique impacts on the trailing hemisphere of the Moon. 

5. The ratio of the number of collisions with the Earth and those with the Moon is found to be 
20.4 + 0.1. 
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